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Abstract 



The pion-nucleon amplitude is considered in the vicinity of the elastic 
scattering threshold within a relativistic dynamical model dressing the uNN 
and ttNA vertices self-consistently with an infinite number of meson loops. 
The dressing is formulated as solution of a system of coupled integral equa- 
tions incorporating unitarity, crossing symmetry and analyticity constraints. 
The calculated scattering lengths and the sigma-term agree with recent data 
analyses. The dressing is important in this model both below and at thresh- 
old. The contribution of the A resonance is discussed, including effects of the 
consistent dressing of the vrA^A vertex. A comparison with the approaches of 
chiral perturbation theory and the Bethe-Salpeter equation is outlined. 
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I. INTRODUCTION 



The pion-nucleon scattering amplitude near the physical threshold is an interesting object 
to study for a number of reasons. At the threshold point itself, the amplitude is proportional 
to the s-wave scattering lengths, whose values are known to be strongly constrained by chiral 
symmetry [0. For the amplitude below threshold, one can establish other chiral low-energy 
theorems involving such quantities as the nucleon sigma-term and thus related to the 
pattern of the explicit chiral symmetry breaking of QCD ^ . To extract the sigma-term from 
scattering data one usually analyses the amplitude at the Cheng-Dashen point. Although 
corresponding to unphysical kinematics, this point is of special importance because both 
pions are on-shell there and the difference between the amplitude and the sigma-term is 
minimal 

Chiral perturbation theory has been used to study loop corrections to the low-energy the- 
orems and, in particular, to calculate the sigma-term However, the near-threshold region 
is significantly affected by the presence of singularities which may make a non-relativistic 
perturbative expansion unreliable 0. In general, the pion-nucleon amplitude is ill-defined 
(and hence non-analytic) at the Cheng-Dashen point and threshold is a branch point dictated 
by unitarity. To obey unitarity for pion-nucleon scattering exactly, one can solve the rela- 
tivistic Bethe-Salpeter equation with a tree-level potential; this yields a good description of 
the phase shifts and scattering lengths as well as allows one to calculate the sigma-term [jlO 



At the same time, the models based on the Bethe-Salpeter equation usually do not preserve 
crossing symmetry (however, see [|ll|) which plays an important role in the derivation of the 
low-energy theorems P]. 

In this paper the pion-nucleon amplitude is studied in the near-threshold region using a 
relativistic dynamical model which incorporates essential constraints from unitarity, analyt- 
icity and crossing symmetry. The effective lagrangian of the model includes pions, nucleons, 
the A resonance, the p and a mesons. The main distinguishing feature of this approach is a 
special method of calculating pion-nucleon and other meson-baryon loop corrections to free 
propagators and bare vertices. An infinite series of loops is summed up by solving a system 
of coupled integral equations for the dressed vertices and propagators. These equations are 
formulated so that constraints from unitarity, crossing symmetry and analyticity are ful- 
filled. This formalism was expounded in Refs. with a simplified treatment of the A 
resonance. In the present paper, not only the irNN vertex and the nucleon self-energy, but 
also the ttA^A vertex and the A self-energy are dressed up to infinite order. The dressing 



procedure of |]T2[ was extended in |T^ to include photons while preserving gauge-invariance. 
Unitarity of the coupled-channel S-matrix above threshold was ensured since the dressing is 
consistent with the application of the dressed vertices and propagators in the K-matrix ap- 
proach. In this way a good description of intermediate-energy pion-nucleon scattering, pion 
photoproduction and Compton scattering was obtained, and the nucleon electromagnetic 
polarisabilities were evaluated and found to be in agreement with experiment |T^JT5|| . Since 
all parameters of the model (resonance coupling constants and a regularising cutoff) were 
fixed at the intermediate energies in Ref. |T^, the present calculation of the pion-nucleon 
amplitude below and at threshold is determined solely by the loop dynamics. 

The outline of the paper is as follows. The pion-nucleon amplitude at threshold and at 
the Cheng-Dashen point is defined in Section where also relevant low-energy theorems are 
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cited. The formalism of the model is described in Section |T|, including the integral equations 
for the dressed vertices and propagators and the method of their solution. In Section 
the invariant pion-nucleon amplitudes in the near-threshold region are explicitly expressed 
via the dressed functions. The main results of the paper will be presented in Section 0. In 
particular, we will examine effects of multiple meson loops, of the A resonance and of the 
a and p mesons, both below and at threshold. The role of the A will be discussed further, 
including the effects of the dressing of the vriVA vertex. We will argue that the analyticity 
constraints incorporated in the dressing procedure are essential for the description of both 
the scattering lengths and the subthreshold coefficients in the same dynamical approach. 
Our formalism and results will be compared with the approaches of chiral perturbations 
theory and the Bethe-Salpeter equation in Section Concluding remarks are made in 
Section |VI]. 

II. PION-NUCLEON AMPLITUDE NEAR THRESHOLD 

The standard isospin decomposition of the pion-nucleon amplitude is |[TH] 



M^p = 5^pM+ + i[r„, r^]M-, (1) 
where Tq, are Pauli matrices for the pion isospins. The spin structure of the amplitude is 

M± = u{p') |a± + 1±Ib^~^ uip) = u{p') {Z}± - (2) 

where k and k' {p and p') are the four- momenta of the initial and final pions (nucleons), 
respectively, u{p) is the Dirac four-spinor. The invariant amplitudes A^, B"^, D"^ depend 
on the Mandelstam variables s = {p + k)"^, u = {p — k'Y and t = {k — k')"^. We will also 
use the standard kinematic variables u = {s — u)/{Am) and vb = {t — 2/i^)/(4m), where 
m = 0.939 GeV and fi = 0.138 GeV are the nucleon and pion masses. The two sets of 
invariant functions in Eq. (H) are related a.s D = A + vB. 

In the following we shall calculate the value of the pion-nucleon amplitude at the sub- 
threshold Cheng-Dashen point, i.e. at z/ = 0, t = with both pions on shell, /c^ = /c'^ = /x^. 
One quantity of special interest is the sigma-term, which is related to the pion-nucleon am- 
plitude at the Cheng-Dashen point as 

S = F^2limD+(z/,t = 2/x2), (3) 

where F^^ = 92.4 MeV is the pion decay constant. The bar indicates that the tree-level 
amplitude evaluated with the pseudovector ttNN vertex (usually called the "pseudovector 
Born contribution") is subtracted, 

B^ = B--'-^^^, (4) 
5- = B--»!(^ 1). (5) 



3 



We take the value g = 13.02 |17| for the ttNN couphng constant. The physical masses 
and coupling constants of the particles used in this calculation are summarised in Table |. 
According to a chiral low-energy theorem the sigma-term Eq. (|^) at the Cheng-Dashen 
point equals the scalar form factor of the nucleon up to corrections of order The scalar 

form factor can be related to the explicit chiral symmetry breaking (see, e.g., Refs. @J^). 
Another low-energy theorem MM concerns the amplitude D~ , requiring that the coefficient 



C = 2F' lim 



D {u,t = 2/i2 



(8) 



should approach unity in the chiral limit (i.e. for a vanishing pion mass), up to corrections 
of order 0(/i^). 

In addition to the pion-nucleon amplitude in the subthreshold region, we shall also cal- 
culate the s-wave scattering lengths which characterise the amplitude at the threshold point 

Th = {s = {m + u={m- ^if, t = 0}: 



,1/2 



D+ + 2D- 
47r(l + fi/m) 



(9) 



Th 



,3/2 



47r(l + /i/m) 



(10) 



Th 



corresponding to the total isospins 1/2 and 3/2, respectively. The low-energy theorem |T| 
asserts that at lowest order the numerators in Eqs. (|) and (0) equal /i/ (F^) and —fi/{2F^), 
respectively. 

III. DESCRIPTION OF THE MODEL 

Our calculation of the near-threshold pion-nucleon amplitude is based on the approach 
of Refs. 



12-14 



However, the treatment of the A resonance is significantly improved in 
the present version of the model, as will be explained in more detail below. In this section 
we describe our approach, focusing on the ingredients which are most relevant in the near- 
threshold region. 



Structure of the amplitude 



The ttN amplitude below and at threshold is purely real. It this model it is constructed 
as the sum of the s-, u-channel nucleon and A exchange graphs, plus the t-channel p and a 
meson exchange graphs. 



M = Ms + Mu + M^ + + Mr, 



(11) 



as shown in Fig. |I]. These are not simple tree diagrams, but rather skeleton diagrams as 
they comprise dressed vertices and propagators. Being a solution of a system of coupled 
integral equations, the nucleon and A propagators and vertices are dressed with meson 
loops up to infinite order, while the p and a propagators are calculated in a one vrvr loop 
approximation. Thus the central element of the approach is the calculation of the dressed 
vertices and propagators, which will be described in the following. 



Structure of the dressed vertices and propagators of the nucleon and A 

The irNN vertex required throughout the dressing procedure has only one of the nucleons 
off the mass shell with the other nucleon and the pion being on-shell (the so-called half-off- 
shell vertex). For an incoming off-shell nucleon with the four- momentum squared p^, the 



most general Lorentz- and CPT-covariant structure of such a vertex is [1 



r.m = T^r'lGpsip') + ^p^Gpvip')], (12) 
L 2m -I 

where Gps^pvip^) are pseudovector and pseudoscalar form factors, to be computed below. 
The Lorentz-invariant expression for the nucleon self-energy is written in terms of two in- 
variant functions Aip^) and B{p^): 

T.{p) = ^L{p)-{Z2-l){i>-m)-Z2 5m, ^l{p) = A{p')i) + B{p')m , (13) 

where S/, denotes the loop contributions to the self-energy. The complete self-energy S(p) 
contains also the counter-term contribution with renormalisation constants Z2 and 5m ad- 



justed to provide the correct pole properties Eq. (|32D of the dressed nucleon propagator ||2T 



S{p) = 1 = ^ + ^(P') (14) 

where for later use we have introduced the self-energy functions 

/ 9n ^ , / 9n ^/ 9n mBirP') + Zoim — 5m) 
aip^) = Z2-Aip^), ap') = t 2 -■ 15 

The nucleon self-energy will be computed consistently with the uNN vertex. 
We choose the following form of the ttA^A vertex: 



We use the fully relativistic formalism, with the metric tensor, 7 matrices and other general 
conventions of pO[| . 
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where p and k are the 4-momenta of an incoming A and of an outgoing pion, respectively, and 
mA is the mass of the A. The part Ga(p^) of the form factor in Eq. (|16|) will be calculated 
in the dressing procedure. The real function F^^jy^^ij) — k)^) depending on the nucleon 
momentum is needed for convergence of the procedure. The isospin transition operators 



are defined by the relations |22 



T.Tj = 5„,-^, TiT^ = l. (17) 

To keep the calculations tractable, the vrA^A vertex in Eq. (0) is not chosen in the most 
general Lorentz-covariant form (for comparison, throughout the calculations we maintain 
the most general structure of the tcNN vertex). It is important however that the vertex 



Eq. (16) has the property of "gauge-invariance" [p3[] , 

P-V{k,p) = 0, (18) 

which allows us to eliminate the background spin 1/2 component of the A propagator ||2 
and to keep only the spin 3/2 component 



Sl'ip) - S^(p)^3;2(P) - ^(p2^[p2 J^V)] ^372(P)' (19) 



where the spin 3/2 projection operator 

Formulae completely analogous to Eqs. (|13|,|15D hold for the A self-energy. Although treating 
the A as a pure spin 3/2 state does not improve the description of pion-nucleon scattering 
phase shifts as compared to the conventional treatment |]T0[, it significantly simplifies the 



dynamical calculation of the A self-energy: we need to compute only two self-energy func- 
tions A^Ip"^) and B^lp"^), instead of 10 invariant functions which would be required if 
the spin 1/2 background were not eliminated. 



Integral equations for the dressing and their solution 

The TiNN form factors Gpv,ps{p^)i the nucleon self-energy functions A{p'^),B{p^), the 
ttA^A form factor Ga(p^) and the A self-energy functions A/^[p'^),B/^{p'^) are calculated 
by solving a system of coupled integral equations. This amounts to dressing these two- 
and three-point Green's functions with meson loops up to infinite order. In the earlier 



version of the model |12|-[l^ the A resonance was not treated completely consistently with 
the nucleon: the A self-energy was computed up to one vrA^ loop only and the dressing of 
the nNA vertex was not included. However, considering nucleon Compton scattering, we 
showed that such simplified A dressing, while being generally adequate, can lead to 
problems at low energies. Therefore, in the present work we refine the dressing procedure 
so that the nucleon and A are now treated on the same footing. 

The dressing equations will be formulated using the following notation. A generic Green's 
function Q{q) is a sum of independent Lorentz-structures (e. g. 1, ^, 7^, etc.), each of which 
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is multiplied with a Lorentz-invariant function depending on (such as form factors or 
self-energy functions) . If we use only imaginary or only real parts of the invariant functions, 
the result will be denoted as Qi{q) or Quil), respectively. If Q{q) is calculated from a loop 



integral, then according to Cutkosky rules Qi{(l) is proportional to the discontinuity of 
the integral across the unitary cut in the complex plane (due to pinching poles of the 
propagators in the integrand) and Qniq) is the principal- value part of the integral.^ 

In our case, the principal-value and pole parts of the dressed ttNN vertex are denoted 
as Tji{p) and Tj{p), respectively. The expression for Tji{p) or Tj{p) is obtained by us- 
ing, respectively, only the real or only the imaginary parts of the form factors Gpv,ps{p'^) 
in the right-hand side of Eq. ([T^). The same applies to the ttNA vertex Eq. (|1^): to 
obtain (V^)ij(A;,p) or (y^)i{k,p) we use only ReGA(p^) or ImGA(p^), respectively. Simi- 
larly, the pole part S/(p) of the nucleon self-energy Eq. (|T3|) contains only ImAlp"^) and 
ImS(p2), and the principal-value part S_r(p) only ReA^p"^) and ReB{p'^). The pion propaga- 
tor D(k) = [k'^ — + iO] ^ does not get dressed, therefore its imaginary part comes from the 
on-shell pions: Dj{k) = 6{k^ — n'^)6{ko). In the same way, we retain only the dominant pole 
contribution to the discontinuity of the nucleon propagator: Si{p) = {'p + m)5ijP' —rn?)6{pQ). 
The resonance propagators do not have poles on the physical Riemann sheet, so their dis- 
continuous parts come solely from their self-energies. For example, the discontinuity of 
the dressed A propagator Eq. ([T9|) is obtained by keeping only the imaginary parts of its 
invariant functions: 

( A j/^ ri{p^)[p'^ — ul)'^{p'^)]^ ?7(p2)[p2 — u;2(p2) ] } ^^^^ ' ^ 

and analogously for the p and a mesons (unlike the A, however, the propagators of the 
meson resonances are dressed in a one vrvr loop approximation only, as will be discussed in 
more detail below and in Appendix A). 

With the introduced notation, the system of dressing equations can be written 

r/(p) = ^ fd^k Tr{p' - k)S{p' - k)VR{p' - k)Si{p - k)Di{k)TR{p) 

+ -^Jd'k iV,)ni-k,p' - k)S^^{p' - k){V,)R{-q,p' - k) 

xSjip-k)Dj{k)rR{p) 
-^Jd'k Tr{p' - k)S{p' - k){V,)R{q,p~ k){S^^)j{p - k) 

xDj{k)iV,)R{k,p-k) + rf(p), (22) 



Re r ^ = ro (P ) + - I dp'' ,2 \ ' (23) 

GpS J [ Gps J TT J(m+M)2 p'^ - p' 



^If the theory obeys analyticity (causality) constraints, the pole and principal-value parts of a 



loop must be related to each other through a dispersion integral [27|. 
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Mp) = J d'kSiip - k)Dj{k)TM , (24) 

Re|^|(/) = -/ dp'' J ^ , (25) 

{V,)i{q,p) = ^Jd'k Tnip' - k)S{p' - k)Tn{p' - k)Sj{p - k) 
xDiik){V,)R{k,p) 

+ ^Jd'k iV,)Ri-k,p' - k)Stip' - k)iVx)Ri-q,p' - k) 

xSj{p - k)Dj{k){V,)R{k,p) + {Vpj{q,p) , (26) 



ReG^{p') = Gl{p') + L r ^^.2 imGMP J ^^^^ 

71 J{rn+n)^ p'^ — p^ 

V'f (p) r — 

Sf(p) = -j^ J d''k{V,)n{k,p)Sj{p-k)Dj{k){K)n{k,p) , (28) 
^^U ) (^') = - / .^P" ^ 2 ' (29) 

[BaJ TT J{m+fj.)2 p'^ -p^ 

where V denotes taking the principal-value of an integral, and isospin factors have been 
absorbed in the coefficients on the right-hand side. The inhomogeneities rj'^(p) in Eq. (p2D 
and (V^°^)/('?,p) in Eq. (p6|) contain p and a mesons, as described in Appendix A. The 
real functions G^py psiP^) ^^"^ bare nNN and vrA^A form factors, respectively. To 

see that the system of dressing equations is neither under- nor over-determined, note that 
Eqs. (|22| , p^ and ( p8|) have two independent spinor structures each. Thus Eqs. (|2^-|29D 



are 14 scalar equations for 14 scalar functions ImGpvip'^), Ini(jp5(p^), lmA{p''), Imi?(p^), 
ImGA(p'), ImAA(/), Im5A(p'), ReG'py(p'), ReGpvip'), ReA{p'), ReB{p'), ReG'A(p'), 
Rev4A(/), Re5A(/). 

Formally, Eqs. (p^-p9D constitute a coupled system of nonlinear integral equations. De- 
spite its quite complicated analytic form, this system of equations has a rather transparent 
meaning (see Fig. 0). The equations are solved by iteration, starting with input bare form 
factors G^pypg, G^. In the course of iteration one effectively sums up an infinite series of 
meson-loop corrections to the bare vertices and free propagators. At each iteration step 
we first calculate the discontinuities of the loop integrals through the Cutkosky rules; these 
pole parts are then used in dispersion integrals to compute the corresponding principal- 
value parts. The details of the computation technique can be found in |ll2|jr3|| . Here we will 
recapitulate only the most important points and discuss the new issues arising due to the 
consistent incorporation of the A resonance in the dressing procedure. 

The use of bare form factors G^py pgijP') and is necessary to regularise the equations.^ 



^An attempt to get rid of the bare form factors by using subtracted dispersion integrals in 
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We choose the purely pseudovector structure for the bare vrA^A^ vertex, i. e. G^gijP') = 0, 
since the derivative couphng of pions at low energies is dictated by chiral symmetry. The 
bare nNN from factor is chosen in the form 



G'pvip') = f exp 



ln2 



(30) 



where is a half- width. The bare coupling constant / = /^nn is adjusted so that the 
dressed ttNN vertex is normalised on-shell to the physical ttNN coupling constant: 

limr(p) = ^7 7', (31) 

where the "sandwich" between the spinors of the initial and final nucleons is implicit. The 



usual renormalisation |21| of the dressed nucleon propagator 



lim S(p) = (32) 

]^^m f — m 

is imposed by adjusting the field renormalisation constant Z2 = and the mass shift 
8m = 6mfyf (see Eqs. ([T3|) and (p!4D ). The renormalisation of the TriVA vertex and of the A 
propagator is done similarly, except that now the pole properties are required of a propagator 



with only the real parts of the self-energy functions |26,12|. There are three corresponding 
renormalisation constants: f-^-NA, and 6m Note that due to the coupled nature of the 
dressing Eqs. (^^- p9|) the six renormalisation conditions for the vertices and propagators 
can be obeyed only simultaneously. Thus fnNN, Z2 , Sm^, and fnNA, Z^, 5mA are inter- 
dependent. The complete set of renormalisation constants obtained in the calculation are 
given in Table |l|. 

We stress that the half-width is not a completely independent parameter: the itera- 
tion procedure converges only for A^ < (A^'*'^)^, and it is important that (A^*^"^')^ is much 
larger than the energy scale due to the explicitly included particles. With the set of param- 
eters used in this calculation, (A^*^^)^ ^ 3 GeV^. When a convergent solution does exist, 
it is reached in practice after about 30 iteration steps for A^ ^ (A^"^)^. The bare vriVA 
form factor G%{p\) as well as the form factors F^na{p%), Fp^^{pp), F„^^{pl), FpNN{p%) 
and Ff^NNipj^), appearing in the vertices with the resonances, have the same exponential 
form as Eq. (0), but peak at the masses of the corresponding particles (see Appendix A). 
We assume the width A^ of these bare form factors to be the same for all resonances, and 
set it close to the maximal value allowed by the convergence requirement. This is done in 
keeping with the general emphasis of our approach on the loop dynamics as determined by 
the dressing rather than on fitting additional parameters. Ideally, values of A|, for differ- 
ent meson resonances should come from a dynamical dressing of these mesons on the same 
footing with the nucleon and the A. Such an extension of the model is certainly feasible, 
although has not been done yet. 



Eqs. (^,25,27,29) fails since each new iteration step would require more subtractions than the 
previous. 
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The effects of the loop corrections on the irNN vertex are similar to those discussed in 
detail in |12;14], where the A was not dressed consistently. We mention the main points here. 
The dressing generates an energy-dependent admixture of the pseudoscalar coupling, which 
at low energies remains much smaller than the pseudovector component and becomes more 
prominent only at intermediate energies. It is important that the dressing does not allow for 
a large pseudoscalar admixture to develop in the low-energy region. The pseudovector form 
factor in narrowed by the dressing. This softening persists independently of the functional 
form of the bare form factor (provided the latter falls sufficiently fast at infinity) and is 
stronger for wider bare form factors. 

Like the half- widths and A|j, the coupling constants QpNN, i^p, daNN, dann and /o-tt^ of 
the p and a mesons are mutually constrained by the requirement that a convergent solution 
of the dressing equations should exist. Fig. | shows the area of convergence in the space 
of the coupling constants gaNN and gpNN, the other parameters being fixed at their values 
given in Table |IT1|. The convergence area has a nontrivial shape for small QaNN and large 
gpNN (or for small QpNN and large g^NN)- For these values there is less cancellation between 
loops with cr's and those with p's, which precludes convergence. Generalising the illustration 
in Fig. 1^, convergent solutions of Eqs. (p2H29|) can be found only in a certain subspace of the 
space spanned by A^, A|., gpNN, i^p, gaNN, gu-K-K and fa-K-K- The remaining freedom in the 
parameters was removed by calculating the pion-nucleon phase shifts using the dressed K- 
matrix approach [T^ and comparing them with data analyses at intermediate energies. 
The dot in Fig. ^ corresponds to the thus optimised coupling constants gpNN and gaNN- The 
full set of relevant parameters from are reproduced in Table |T|. We emphasise that 
in fixing these values we did not use any results of data analyses for the scattering lengths 
or for the subthreshold amplitudes. Therefore there are no free parameters in the present 
calculation of the near-threshold amplitudes. 



IV. SCATTERING AMPLITUDE IN TERMS OF THE DRESSED VERTICES 

AND PROPAGATORS 

Having solved the dressing equations, we proceed to evaluate the invariant amplitudes 
and B"^ defined in Eq. (|^). First we write the s- and u-channel diagrams from Eq. (pH]) in 
terms of the dressed nNN, nNA vertices Eqs. (|T2|,|T6|) and dressed nucleon, A propagators 
Eqs. (14|,19). Then we add the t-channel diagrams written in terms of the relevant meson 



vertices and propagators Eqs. ( |X3| jA5|JA 1 0| , |A 1 4| , [A 1 5| , |A2 1| ) . At the considered kinematics, 
only the real parts of the invariant functions of the two- and three-point Green's functions 
enter in the amplitude. After doing some straightforward algebra, the contribution of the 
nucleon, p and a exchange diagrams to the invariant amplitudes can be written 



A' 



+ 



Ghmis) -m] + ^-s{s)GMs) _ ^ G^^^ ^ ^^^^^^^ _ 

ais)[s-e{s)] 

GlsiuMu) + Gps{u)Gpy{u) ^^ _ ^ 9^[m + e(«)][« - 

a{u)[u — ^^(s)] 
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gaNNFanTrit) 



QcnvTrf^ fa 



+ - t 



z'^[t-cm 



(33) 



ais)[s-eis)] 



Gls{umu)-m] + 



Gps{u)Gpv{u) 



m 



[u — w?] + 



Glyiu] 



Aw? 



[m + ^(■u)]['U — w? 



a{u)[u — .^^(s)] 



+ 



gpNNgpTnrl^pFpmr{t){u — s) 

2mZP [t - A2(t)] 



(34) 



o+ m 4m"' 



+ 



G'p5(n)G'py(M) 



m 



B- 



a(u)[u — ^^{u)] ' 



(35) 



m 



Gps{u)Gpv{u) 



m 



[m + ^{u)] + + « + 2me(M)] 

4m^ 



a{u)[u — ^^(u)] 



+ 



(1 + /Cp)F^,,(t) 



ZP[t-A2(t)] 

The A exchange diagrams from Eq. ( pA]) are given by (restoring the isospin indices) 

Gl{s) 



(36) 



Ial3 — I Oq/3 



ri{s) [s — i-o'^(s)] 



^(p) {f + f)k'-{p' + k').k'^. 



(37) 



X [jf^ + ^ + a;(s)] P372(p + k) [{^ + ^)k, -{p + k)-k 7,] u{p) 
(M^\ - (s '^^'^A ^a(^) 

X - r + ^(^)] ^372(P - \-{^ - V)k'. + {P- k') ■ A;' 7.] u{p) . (38) 



Mp') [-{f-m,+ ip'-k)-k^,] 



For brevity we will not give the explicit decomposition of Eqs. ( P7| , PSD in terms of the 
invariant amplitudes. 



Properties of analyticity and crossing symmetry 

Due to the use of the cutting rules and dispersion integrals in the formulation of Eqs. (p^ 
291) , the two- and three-point Green's functions obtained by solving these equations possess 
the correct analyticity structure associated with the nucleon and A exchanges and obey the 
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two-body (vrA^) unitarity. Furthermore, the vrA^ amphtude obeys the crossing symmetry 
requirements that A'^, A~ / {s — u), / {s — u) and B~ be invariant under the replacement 
s u. The crossing is respected due to our using the dressed two- and three-point Green's 
functions in both s- and u-type diagrams. At the same time, the t-channel analyticity 
structure is not fully reproduced because the t-channel cuts are taken into account only 
through the tttt loops in the p and a propagators but the loop corrections to the vertices 
with the p and a mesons are not included. Also, the four-point one-particle irreducible 
diagrams-such as the box graph-are not included in the dressing. 

The omitted dressed one-particle irreducible four-point diagrams can be thought of as 
being of order O{o?) in a certain formal expansion [Q, where the parameter a characterises 
the level of analyticity violation in the model. In this expansion, the lowest order 0{oP) 
corresponds to an amplitude with no dressing, in which case the violation of analyticity 
in maximal. The next order 0{a}) corresponds to an amplitude in which the one-particle 
reducible (with respect to the s-channel cuts) graphs contain the dressed propagators and 
vertices. Thus, at order 0{a}) analyticity is restored at the level of two- and three-point 
Green's functions, as is done in present dressing procedure. The higher orders in a are 
described by induction in terms of reducibility of n-point Green's functions constituting the 
amplitude. 

The parameter a can be defined as follows. If a scattering amplitude T(a;) of a process 
can be represented at small energies a; as a power series 

T{ijj) = Co + ciuj + C2uj'^ + . . . , (39) 

then the coefficients Cj could, in principle, be computed using the same model but in two 
different ways: 

1. Low-energy {LE) evaluation: compute Ci{LE) by evaluating T{uj) directly at low 
energies; 

2. Sum-rule (SR) evaluation: compute Ci{SR) by calculating appropriate total cross sec- 
tions and integrating them in the sum rules corresponding to each particular coefficient. 

These two ways of evaluation should give identical coefficients Ci{LE) = Ci{SR) = Ci pro- 
vided the analyticity of the model is exact. In practical calculations there will always be 
discrepancies between the two methods, which can be used to estimate the violation of 
analyticity in the model. Therefore one can quantify the parameter a as 

a^\ciiLE)-CiiSR)\, (40) 

where the choice of a particular coefficient from the series Eq. (|39D could be decided by 



additional considerations. This idea was tested for nucleon Compton scattering in Ref. [[15 
where the coefficients Ci and 02,3 were related to the anomalous magnetic moment and to 
nucleon polarisabilities, respectively. The proposed formal expansion in the parameter a 
can offer a systematic way of improving analyticity properties of dynamical approaches 
applicable at low and intermediate energies. 



12 



The properties of analyticity and crossing symmetry are crucial for the model to provide 
a good description of the amplitude both below and at the physical threshold.^ To study 
this in more detail, we shall focus below on the sigma-term, the Adler-Weisberger coefficient 
C, as defined by Eqs. (|,|]), and on the scattering lengths, as defined by Eqs. (|9|,p!0D. We 
wiU collectively call these quantities the "near-threshold coefficients" . 



V. NEAR-THRESHOLD COEFFICIENTS 



The sigma-term and coefficient C at the Cheng-Dashen point 



On expanding the amplitudes in Eqs. (|33|-|38D around the Cheng-Dashen point s 



u 



rn^.t = (with k"^ = k'"^ = /i^) and using the definitions Eqs. (|],|D, we obtain exphcit 
formulae for the sigma-term and for the coefficient C in terms of the dressed vertices and 
propagators: 



S = -P 



C 




QaNN g 



+ 



dpNN 5'p7r7r(l 




(41) 



) FpTTTT ( 2/i 



+ Ca, 



(42) 



where Sa and Ca contain the effects of the A exchange in the s- and u-channel diagrams 
(see Appendix B). In Eqs. (^T]) and (^) we have made use of the relations 



e(m2 



aim 



2ma{im? 



9, 



(43) 



which follow from the renormalisation conditions Eqs. ( pT] , p2D . 

Eqs. (PI) and ( ^2]) show that the subthreshold parameters S and C are sensitive to the 
values and derivatives of the dressed irNN form factor Gpsip"^) and of the nucleon self- 
energy function a(p^) at the nucleon pole. Since the pseudoscalar form factor in the dressed 
ttNN vertex is very small in the vicinity of m^, the nucleon contribution in Eqs. (^,^2|) is 
much smaller than that of the a and p mesons. The a exchange in the t-channel gives a 
dominant numerical contribution to the sigma-term due to the presence of the non-derivative 
component ~ gamrf^ in the cxTrTr vertex Eq. ( |A14| ). The p meson plays a similar role for the 
coefficient C. It should be pointed out, however, that dissection of Eqs. (41 
"nucleon contribution" , "a and p meson contributions" and a "A contribution 



mto a 
can only be 

regarded as formal here: as already mentioned, Eqs. (|2^-p9D not only determine the nucleon 
and A dressing, but also strongly constrain the other parameters of the model. For example. 



similar conclusion was reached in the framework of the relativistic baryon chiral perturbation 
theory |p, where it was pointed out that the standard low-energy expansion does not reproduce 
the correct analyticity structure in the vicinity of singularities. 



13 



an arbitrary change of the values of the p and a couphng constants would formally change 
the dominant contributions in Eqs. (^ , |42|) , but with these altered coupling constants the 
dressing procedure might not converge at all! 

In what follows we will discuss various ingredients of the dressing by comparing the 
"Dressed" and "Bare" calculations. The former contains the full dressing with the meson 
loops whereas in the latter the bare vertices and free propagators have been used. Note that 
since the pTrvr, pNN, airn and aNN vertices do not get dressed in the model, in both calcula- 
tions they are equipped with the bare form factors -FttataIPat), Fp^Tvipf,), F^mripl), Fpnn{Pn) 
and F^N^{pj^), as defined in Appendix A. The obtained values of the near-threshold coeffi- 
cients are summarised in Table Results of several data analyses are quoted in the last 
row. Various ingredients of the fully dressed calculation are given in the other rows and will 
be discussed below in more detail. 



Pion-nucleon scattering lengths 

The scattering lengths are evaluated by substituting the explicit expressions for the 
amplitudes Eqs. (p^-|38|) into Eqs. (H) and (|1^). The obtained values are listed in Table 
[rV| for the different calculations considered. The large effect of the dressing on a^/^, in 
comparison with the small effect on a^/^, is a reflection of the nucleon s-channel graph being 
influenced by the dressing more than the u-channel graph. 



Dependence on the bare form factor 

To see how our results depend on the choice of the bare form factor Gpy in Eq. (pOD, 
we did a calculation in which the width of the bare form factor was set to A^ = 2.8 
GeV^ (i. e. near the upper limit (A^)™"^ dictated by the convergence) and all the other 
parameters were kept as given in Table |TTJ As was shown in [T^, despite using a much 
wider form factor, such a calculation leads to intermediate-energy phase shifts which are 
similar to those obtained in the basic calculation with A^ = 1.8 GeV^. With this variation 
of the bare width, the sigma-term and the coefficient C change by less than 3%. The 
sensitivity of the scattering lengths is similarly small. The important point here is that even 
using quite different bare form factors the model yields results which are comparable with 
data analyses both at the Cheng-Dashen point and at threshold. As long as the iteration 
procedure converges, the loop dynamics depend only weakly on the details of the bare form 
factor. The usage of the exponential bare form factor is also not essential and dipole-like 
bare form factors lead to similar dressed vertices. By contrast, the role of the dressing is 
quite significant. This can be seen by comparing the rows labelled "Dressed" and "Bare" in 
Table ITVl. 



Effects of the dressing of the a and p propagators 

In the fully dressed calculation, the self-energies of the a and p mesons are dressed with 
one Tin loop, as detailed in Appendix A. By using the free propagators instead, we obtain 
the results listed in the rows "Free a" and "Free p" in Table ^V\. The p exchange contributes 
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solely to the isospin-odd amplitudes while the a exchange acts in the isospin-even channel, 
which is explicitly shown in Eqs. (13^^361). Hence the sigma-term E (being calculated from 
in Eq. (^) is oblivious to the treatment of the p meson while the coefficient C (being 
calculated from D~ in Eq. (|)) is not affected by the a meson. 



Role of the A resonance 

A pole contribution 

By the A contribution to the pion-nucleon amplitudes one usually means the contri- 
bution of the s- and u-channel A exchange diagrams. This definition may be used for a 
comparison of our results with those of the chiral perturbation theory, where the A is usu- 
ally not included as an explicit field in the lagrangian and thus does not appear in the 
loops p2| , |33| ,P| -PI By comparing the rows "Dressed" and "No A poles" in Table we see 
that the pole contribution of the A is small both below and at threshold (the explicit formu- 
lae are given in Appendix B). For example, it changes the sigma-term by 0.26 MeV, which is 
somewhat smaller than the values typically found in other calculations [[^,p|,p5[. However, a 



quantitative comparison of our results with other approaches should be carried out carefully 
since typically one retains a spin 1/2 background in the A propagator (see, e. g., P^J5[]). 
By contrast, in this model we deal with the pure spin 3/2 A due to the gauge-invariant 
structure of the ttNA vertex Eq. 



Effects of the dressing of the irNA vertex 

In addition to the s- and u-channel exchanges, the A resonance enters in the loop correc- 
tions to the ttNN and nNA vertices both of which are dressed up to infinite order. Such a 
contribution of the A through dressing has not been considered before in the context of the 
near-threshold vrA^ amplitude. Due to the coupled nature of Eqs. ([2^-P^, the A dressing 



affects the ttN amplitude directly as well as through its effects on the ttNN vertex and 
nucleon propagator. 

In the subthreshold region, the value of the sigma-term is slightly decreased by the A 
dressing, as can be seen by comparing the rows "Dressed" and "Bare A" in Table (in the 
latter calculation the bare vrA^A vertex and the free A propagator were used). Comparison 
of these results with the column "No A poles" shows that the dressing moderates even 
further the small pole contribution of the A to the sigma-term. Notably, the effect of the 
consistent dressing of the vrA^A vertex is to decrease the sigma-term and to increase the 
coefficient C, whereas the A exchange has the opposite effect. This suggests that meson 
loops should be treated with special care in dynamical calculations of the A contribution to 
the subthreshold ttN amplitude. 



^There exist approaches [Q| including the A expUcitly in chiral lagrangians, in which case the 
A does enter in the loops. Such extensions of the standard chiral perturbation theory will not be 
discussed here. 
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Effects of multiple loops 



The two- and three-point functions obtained as a solution of Eqs. comprise an 

infinite series of meson-loop corrections. The effect of multiple meson loops on the near- 
threshold amplitude can be seen from the row labelled "One loop" in Table |^ This is a 
calculation in which the ttNN and vrA^A vertices, as well as the nucleon, A, p and a self- 
energies, include only one-loop corrections. These vertices and self-energies are obtained by 
iterating Eqs. (p^-pQ]) only once after discarding from them the effective two-loop contribu- 
tions (the latter are the integrals containing one of the cut resonance propagators {S'^)!, 
{D^^)i or {D„)j). To focus on the genuine effects of the loops, no parameters were readjusted 
in this calculation, except the renormalisation constants f^NN, fwNA, , Z^,5miy, 5m^. 
Note that in this case the ttN scattering amplitude is not exactly an amplitude including 
only one-loop corrections to the tree-level approximation. Nevertheless, taking into account 
that most of the effects of the dressing are due to the vertices, the difference between this 
"one loop" and the fully dressed calculations can serve as a reasonable estimate of the ef- 
fects beyond few perturbative loop corrections in the amplitude. Table |^ shows that the 
multiple-loop effects are small, being comparable in magnitude with the contribution of the 
A resonance (Ref. gives an upper limit of ~ 2 MeV for the A contribution to the sigma- 
term). The value of the sigma-term comes mostly from the one-loop calculation, whereas the 
multiple loops and the A give contributions which are beyond the precision of the present 
data analyses. 



VI. COMPARISON WITH THE APPROACHES OF CHIRAL PERTURBATION 
THEORY AND THE BETHE-SALPETER EQUATION 

The near-threshold parameters obtained in this model are compared in Table |V| with 
calculations in chiral perturbation theory (xPT) and with those based on the Bethe-Salpeter 
equation (BSE). We quote the results of the heavy-baryon and relativistic (with the infrared 
regularisation) baryon formulations of chiral expansions, referred to as HBxPT and RB^PT, 
respectively. There is some similarity between the present model and the other approaches, 
as well as a number of differences. In principle, it would be desirable to compare the loops 
calculated in this model with the chiral loops and with the loops generated by the BSE at 
the most elementary level, i. e. by analysing the expansions of a simple Green's function 
such as the nucleon self-energy. Such a comparison is hardly meaningful, however, since 
in general the Green's functions are model- and representation-dependent in perturbative 
as well as nonperturbative approaches (see, e.g., ||3^ and references therein). Instead, the 



values of the near-threshold parameters assembled in Table |V| can be used to illustrate the 
most important similarities and differences between the four approaches: the present model, 
HBxPT, RBxPT and the BSE. 

The basic dynamical degrees of freedom in all these approaches are the nucleon and the 
pion. However, unlike the two chiral formalisms, the present model and the BSE include 
the A resonance and the p and a mesons explicitly. In general, these degrees of freedom 
give important contributions to the dressed Green's functions. In chiral approaches effects 
of the resonances are typically encapsulated in low-energy constants using a resonance- 



saturation hypothesis [p4|32| (as mentioned above, we do not touch upon extensions ||34 
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of the chiral expansion with the A as an exphcit degree of freedom). In our model, as 
well as in xPT, the ttNN vertex is predominantly pseudovector at low energies. A small 
pseudoscalar component in the vertex is associated with an explicit symmetry breaking. In 
all four approaches, the nontrivial dependence of the amplitudes on kinematical variables 
is determined by loop expansions. The methods of organising the expansions are different 
however, the methods of evaluating the loop diagrams. 

The chiral expansions are series in small pion mass and momenta [^. In the presence of 
a nucleon, a one-to-one correspondence is established in HBxPT between the power of each 
term in the expansion and the number of pion loops required for its calculation PU|. This 



nonrelativistic expansion is well-behaved in the part of the low-energy region which is not 
too close to the singularities of the amplitude dictated by unitarity and analyticity. However, 
it is necessary to rearrange the chiral expansion in order to calculate the amplitude near the 
singularities. This is most efficiently done in RB^PT by using an infrared regularisation 
in the relativistic formalism and representing the amplitude through dispersion integrals. 
(Within HBxPT, the correct singularity structure can be captured only by summing an 
infinite number of loop corrections.) The main emphasis of the approach of the BSE is on 
obeying unitarity in the relativistic formalism |^ , p!0| . The necessary loops are summed up 
effectively by solving the linear integral equation for the scattering amplitude. However, due 
to the use of a simplified kernel in practical calculations, the BSE does not yield a crossing 
symmetric amplitude. 

The unitarity, crossing and analyticity constraints are used as the governing principles 
for organising the meson-loop expansion in this model. The (nonlinear) dressing equations 
represent a self-consistent procedure of using the cutting rules and dispersion relations. 
By solving them we effectively sum up an infinite series of loop diagrams in such a way 
that the essential singularity structure of the two- and three-point Green's functions is 
correctly reproduced. However, in contrast to the chiral expansions or the BSE, one-particle 
irreducible four-point loop diagrams (such as the triangle or box graphs discussed in [^]) are 
not calculated explicitly in our model, which entails a violation of analyticity structure in 
the t-channel. Also, in this model we need to compute the dressed vertices only in the s- 
and u-type diagrams, i. e. for on-shell external pions only. Consequently, the dependence 
of the amplitude on the momenta of external pions is not obtained. For instance, to obey 
the Adler consistency condition D {u = 0,t = fi"^) = (with either k"^ = /x^, /c'^ = or 



fc^ = 0, fc'^ = yU^), the aNN vertex and the box graph will have to be dressed in this model, 
thus requiring an extension to order O(a^) in the "analyticity violation expansion" outlined 



in Section [V 



In the calculation of RB^PT the bulk of the sigma-term at the Cheng-Dashen point 
comes from a contact four-point vertex. This contact vertex is proportional to the pion 
mass squared and thus explicitly breaks chiral symmetry of the lagrangian. In the present 
approach, the largest contribution to the sigma-term is due to the t-channel a exchange 
with the component ~ gaimfJ' in the airir vertex Eq. ( |A14|) . This non-derivative part of the 
a exchange does not vanish in the chiral limit and is therefore analogous to the symmetry 
breaking contact term of RB^PT. It should not be concluded, however, that the loop con- 
tributions are of minor importance in the near-threshold region: it is through the loops that 
the essential unitarity, analyticity and crossing constraints are incorporated in this model 
as well as in RB^PT. 
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Although the means of regularisation of the loop integrals employed in this model and in 
the BSE (usage of bare form factors) are different to those utilised in the RB^PT (infrared 
and dimensional regularisation), they are similar in that they generate spurious singularities 
of the amplitude. These unwanted singularities should be removed from the region of physical 
interest. For example, the amplitude in the RB^PT has an unphysical pole at s = 0, which 
should be safely far from the relevant near-threshold region Similarly, the bare form 
factor in this model is rather wide, ensuring remoteness of its singularities. There is an 
important difference between the evaluation of the loop integrals in our approach and in the 
BSE. In the latter method the loops are computed using such techniques as, e. g., the Wick's 
rotation and Feynman parametrisation, whereas in our framework the loops are calculated 
through the successive application of the Cutkosky rules and dispersion relations. While 
being equivalent in local field theories, these two methods of loop evaluation are likely to 
differ when the loops are regularised by form factors. Therefore the analytical properties of 
the amplitudes dressed in the BSE are probably different to those generated by our approach, 
which may have important consequences in the near-threshold region. 

VII. CONCLUSIONS 

A consistent dynamical calculation of the pion-nucleon amplitude in the near-threshold 
region should serve as a bridge between the physically accessible low-energy data and chiral 
low-energy theorems reflecting the QCD dynamics in the nonperturbative regime. In this 
paper we have shown that essential analyticity constraints can be incorporated in a self- 
consistent dressing procedure, resulting in a reliable description of the amplitude at the 
Cheng-Dashen point and at threshold. In particular, the pion-nucleon scattering lengths, 
the nucleon sigma-term and the Adler-Weisberger coefficient C evaluated in our approach 
are all consistent with the recent data analyses. 

The values of the near-threshold coefficients depend crucially on the treatment of the 
a and p meson exchanges. The approach also includes a consistent dressing of the vrA^A 
vertex and of the A propagator. This allows us to study the role of the A for the nucleon 
sigma-term in more detail than has been done previously. In particular, we have found that 
the contribution of the A in the near-threshold region should be considered on a par with 
the effects of multiple loops. 

This dynamical model suggests that effective approaches incorporating the constraints 
of relativistic invariance, unitarity, crossing symmetry and analyticity can reveal important 
aspects of the low-energy strong interaction. 
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APPENDIX A: CONTRIBUTIONS OF THE p AND a TO THE DRESSING 

EQUATIONS 



The term T''j"{j)) on the right-hand side of Eq. (p2D comprises 6 loop diagrams including 
the p and a degrees of freedom^, 



rr(p) = E { (r?)n(p) + (rj)n(p) }, (ai) 

n=l 

where the different terms are described in the following. 

(r?)i(p) = ^1 d'kTf''{k-q,p-k)Si{p-k)TMD^;{q~k)Tr{-k.q)Di{k), (A2) 
where the pim vertex has the structure 



qy-qv + 



(g + q') 



FpU{q + q)) 



(A3) 



In order that a converging solution of the dressing should exist, the pun vertex is equipped 
with a form factor Fpj^T^{{q + q'Y) whose form is the same as that of the bare nNN form 
factor in Eq. (|30D , 



PpTTTV ( k 



exp 



-ln2 



- mlf 



(A4) 



where the half-width A|. is given in Table |IT1| and discussed in Section |1]. By analogy with 
the treatment of the A resonance (see Eq. (p!8D), we choose a gauge- invariant form of the 
pTTTT vertex and therefore retain only the spin 1 part of the p propagator 



D^;{k) 



ZP[k^-X^{k^ 



(A5) 



with the spin 1 projector 



(A;) = gf"" - 



(A6) 



The p self-energy function A^(fc^) is calculated from one vrvr loop as follows. 



nip - Snip + 



ZP 



(AT) 



where HKp ) is the nn loop contribution to the self-energy and Z^, drrip are renormalisation 
constants adjusted to ensure the correct pole properties of the dressed p propagator Eq. ( |A5D . 
As all loop integrals in the model, this loop is evaluated using a dispersion relation: 



^In this appendix we use the notation introduced in Section LD. 
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1mm {k') = J d'q r^--(g, k - q)Dj{q)Dj{k - g)!^ - k, -q) 



Ren^fc^) = - dk 

IT Jau'^ 

The pNN vertex is chosen as 



4aj2 k'^ — k"^ 



2m 



(A8) 



(A9) 



(AlO) 



where the regularising form factor Fpj^]^{p^) depends on the four- momentum squared of an 
off-shell nucleon and has the form of Eq. ( [A^ ) in which m-Ar is substituted for irtp. Describing 
Eq. (^) further, 



iVMp) = g/ d'kSj{p-k)TP^''''{-k,p)D{k-q)Tr{-qA-k){D;^)i{k), (All) 
where Eq. (0) has been used; 

(r?)3(p) = ^ / rf'A^rf - k)S{p' - k)Tn{p' - k)S,{p - k)Tf'^{-k,p){Dpj{k) . 



(A12) 



The loop integrals with the a meson are similar in structure to those with the p. 

dcrNN 



(ri)i(p) = -^-^Jd'k Sj{p - k)TR{p)DM - k)T''-^{-K q)Dj{k), (A13) 



where the (Ttttt vertex is chosen as 



Qcnrnf^ ~l~ fa 



Fanniiq + q 



f\2\ 



(AM) 



The a propagator D„{k) is obtained from one tttt loop: 



Renj(p2^ 



c'(r) =mi- + 



\mW{e) = I A r^Mg, k - q)Dj{q)Dj{k - g)r— (g - k, -q) 



Ren2(r) 



The two other loops with the a read 



dk 



(A15) 

(A16) 
(A17) 
(A18) 
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(V^Mp) = 0/ d'kSj{p-k)T^''''{-k,p)D{k-q)T^-{-q,q-k){D^)j{k), (A19) 



(rj)3(p) = ^J d'kr'^''''{k,p' - k)S{p' - k)rn{p' - k)Sj{p- k)r'^''''i-k,p){D^)jik), 

(A20) 

where the aNN vertex has the simple structure 

r'^^^(A;,p) = g.NNF^NNip'), (A21) 

with the form factor FrjNNip^) = Fp^i^^p^). 

The term (V^'^)/(q',p) in Eq. ( PB| ) comprises two loop integrals whose form is given by 
Eq. ( |A2|) and Eq. ( |A13|) in which the irNN vertex Tfi{p) is replaced with the ttNA vertex 

{V,)R{k,p). 

APPENDIX B: CONTRIBUTION OF THE A AT THE CHENG-DASHEN POINT 

The contributions of the A resonance to the sigma-term and to the coefficient C in 
Eqs. (^,^) are obtained by decomposing Eqs. ( P7|J5B| ) into the invariant amplitudes Z)^ and 
evaluating the latter at the Cheng-Dashen point. The expressions in terms of the dressed 
ttA^A form factor and the dressed A self-energy functions (as defined in Eqs. (|TB|,|IP|)) are 



_ %F^mu{m^)G\{m') ^ 



+ 2P 



d m\ri{m?)[rn? — uj'^{m?) j 



2 



_^ 2mGi(m^)[m^ + co'^(m^) ] [r7(m^)co''(m^) - cj(m^)?7'(m^) ] I 4 
explicitly showing the suppression by powers of the pion mass /x. 
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TABLES 



TABLE L Particle masses (in GeV) and coupling constants used in this calculation (the same 
values were used in Ref. [Q). The masses and the coupling constants of the A, p and a correspond 
to the experimental widths and positions of the resonances as given in [l^. 



m = rriN 


= 




Ulp 




9 = 9-kNN 




QpTVTT 


0.939 


0.138 


1.232 


0.77 


0.76 


13.02 


19.76 


6.07 



TABLE II. The renormalisation parameters: bare coupling constants, mass shifts (in units of 
GeV) and field renormalisation factors. 



/ = firNN fnNA 


6m = 6m N 6m ^ 


6m^p 


6mi = Zi' 




ZP 


Z" 


10.75 21.75 


-0.075 -0.120 


-0.089 


-0.605 0.80 


1.16 


1.17 


1.14 



TABLE III. Values of the parameters of the model, as fixed by calculating the intermedi- 
ate-energy pion-nucleon phase shifts in the dressed K-matrix approach of Ref. . No parameters 
were readjusted in the present calculation. 







QpNN 




QaNN 




fcTTTTT 


1.8 


1.0 


7.0 


2.3 


34 


1.7 


1.8 
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TABLE IV. The pion-nucleon sigma-term, the Adler-Weisberger coefficient C, evaluated at the 
Cheng-Dashen point from Eqs. (|3|,^), and the s-wave scattering lengths, evaluated from Eqs. (p|,|lO|). 
The rows represent the following calculations. "Dressed": fully dressed calculation; "Bare": bare 
calculation, i. e. using the free propagators and no loop corrections to the bare vertices; "Free 
cr": full calculation, but using the free a propagator in the t-channel exchange; "Free p": full 
calculation, but using the free p propagator in the t-channel; "No A poles": full calculation, but 
without the s- and u-channel A exchange pole diagrams; "Bare A": full calculation, but using the 
bare vrA^A vertex and the free A propagator; "One loop" : calculation in which the nucleon and A 
self-energies as well as the irNN and vrA^A vertices are computed up to one-loop corrections only; 
"Data" : results of various data analyses. 





S (MeV) 


C 






Dressed 


73.99 


1.16 


0.175 


-0.087 


Bare 


127.78 


1.31 


0.204 


-0.088 


Free a 


126.41 


1.16 


0.183 


-0.080 


Free p 


73.99 


1.40 


0.210 


-0.105 


No A poles 


73.73 


1.21 


0.175 


-0.087 


Bare A 


74.05 


1.14 


0.175 


-0.087 


One loop 


71.98 


1.21 


0.181 


-0.093 


Data 


64 ± 8 O 


1.15 ±0.02 |29| 


0.173 ±0.003 [p^ 


-0.101 ±0.004 |2^ 




79 ± 7 [0 


0.175 ±0.004 [31] 


-0.085 ±0.027 |311 




71 ±9 









TABLE V. Comparison of the near-threshold parameters evaluated in the present model with 
results obtained in chiral perturbation theory and in the approach based on the Bethe-Salpeter 
equation. The third and fourth order HB^PT calculations |^] presented fits to three different 
phase-shifts, using Ref. to relate the sigma-term and threshold parameters. The RB^PT calcu- 
lation used the data analyses of [29| as input. The BSE results are from Ref. |10|. 





This model 


HBxPT 


HBxPT C(p^) 


RBxPT 


BSE 


S (MeV) 


73.99 


69 
91 

93 


73 
85 
104 


61 


23.6 


C 


1.16 


1.10 
0.82 
1.09 




1.13 






0.175 


0.171 
0.159 
0.175 


0.171 
0.159 
0.176 


0.175 


0.177 




-0.087 


-0.101 
-0.072 
-0.086 


-0.100 
-0.073 
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FIGURES 




FIG. 1. Graphical representation of the pion-nucleon amphtude, corresponding to Eq. (|TT|). 
The single solid lines are nucleons, the double lines are A's, the dashed, zigzag and dotted lines 
are pions, p's and cr's, respectively. The empty and hatched circles denote the ttNN and vrA^A 
vertices, respectively. The propagators and vertices are dressed with meson loops as described in 
Section |l[ 
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H o/-4=f-^ H 



A 



Disp. 



B 



Disp. 



C 



Disp. 



D 



Disp. 



Disp. 



Disp. 

FIG. 2. Graphical representation of the system of integral equations for the dressed two- and 
three-point Green's functions. The notation for the propagators and vertices is as in Fig. For 
each particle, the dressed and free propagators are denoted by thick and thin lines, respectively. 
The triangles denote the counterterms needed to fulfill the renormalisation conditions such as 
Eq. (^2|). The slashes through the loops and the integral signs indicate the use of the Cutkosky 
(cutting) rules and dispersion integrals in the iterative solution of the equations. The outgoing 
nucleons in the vertices (as well as the pions) are on-shell, which is denoted by the crossed lines. The 
correspondence with analytic equations is A (22, 2^, B ^ (24, 25), C <-> (p^,|27|), D 
E ^ (|3|,g^, F ^ ( [HI , [All) . 



(28,2' 
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A ^pNN 



'aNN 







20 40 60 

FIG. 3. Illustration of the interdependence of the coupling constants g^NN and QpNN (the 



other parameters being fixed as in Table piID due to the requirement of convergence of the dressing 
procedure. The area of convergence is sketched using 50 test solutions of Eqs. (p^-p9[), but ignoring 



vrA^ phase shifts. The dot corresponds to the phase shifts fit obtained in Ref. |14| 
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